Theoretical Analysis of the Characteristic Impedance in 
Metal-Insulator-Metal Plasmonic Transmission Lines 
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We propose a closed form formulation for the impedanee of the Metal-Insulator-Metal (MIM) plasmonic 
transmission lines by solving the Maxwell's equations. We provide approximations for thin and thick insulator 
layers sandwriched between metallic layers. In the case of very thin dielectric layer, the surface waves on both 
interfaces are strongly coupled resulting in an almost linear dependence of the impedance of the plasmonic 
transmission line on the thickness of the insulator layer. On the other hand, for very thick insulator layer, the 
impedance does not vary with the insulator layer thickness due to the weak-coupling/decoupling of the surface 
waves on each metal-insulator interface. We demonstrate the effectiveness of our proposed formulation using 
two test scenarios, namely, almost zero reflection in Tee-junction and reflection from line discontinuity in the 
design of Bragg reflectors, where we compare our formulation against previously published results. © 2013 
Optical Society of America 

OCIS codes: 240.6680, 250.5403. 



Surface Plasmon Polariton (SPP) is a TM polarized 
electromagnetic wave traveling along metal-dielectric in- 
terfaces. It can also be explained as an oscillatory mo- 
tion of electrons bounded to the metal-dielectric inter- 
face, where the corresponding field profiles decay expo- 
nentially into the neighboring media from their maxima 
occurring at the interface [1]. SPP can confine light to the 
lateral dimension in the order of tenth of the wavelength 
overcoming diffraction limit for travelling waves in the 
Metal-Insulator-Metal (MIM) transmission lines [1, 2]. 
MIM lines have been recently used to build microwave- 
inspired plasmonic devices, such as filters, branch-line 
couplers, and resonators [3,4], as well as optical-inspired 
devices, such as Bragg refiectors, plasmonic photonic 
crystals, and Mach-Zehnder couplers [5-7]. An effec- 
tive design of the aforementioned devices requires a re- 
liable and accurate method to derive the impedance 
of the MIM line. To the best of the authors' knowl- 
edge, there are only a couple of methods to charac- 
terize the impedance of the MIM line in the literature 
(cf. [8,9]). However, those methods do not accurately 
explain the characteristic impedance of the MIM struc- 
ture for generally designed MIM lines of various thick- 
nesses. The characterization of the impedance of the 
MIM can be useful for a better understanding of the 
MIM lines. It also provides a design tool to generate dif- 
ferent microwave/optical-based plasmonic structures. 

In the following, we propose a formulation for the 
impedance, which is derived by solving the Maxwell's 
equations. This formulation is valid for any thickness of 
the dielectric layer ranging from very thin to very thick 
layers. Then, two design scenarios arc discussed. A design 
of a Tee-junction for minimum reflection using the pro- 
posed formulation and its comparison to other published 
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Fig. 1. Real (blue solid line) and imaginary (red dashed line) 
parts of the effective refractive index of an MIM structure 
(shown in the inset) as a function of the dielectric layer thick- 
ness. The operating wavelength is assumed to be 1.55fxm. 

formulation methods is presented. The second scenario 
analyzes the minimum reflection from a line discontinu- 
ity using the proposed formulation. This line discontinu- 
ity can be utilized in the design of Bragg reflectors. 

Starting with plane- wave excitation of the MIM struc- 
ture and solving the Maxwell's equations, the only prop- 
agating mode is the transverse magnetic (TM) mode. 
The effective refractive index can also be calculated using 
the solution of the Maxwell's equations [I]. The effective 
refractive index for a silicon layer sandwiched between 
silver layers is shown in Fig. I. The real and imaginary 
parts arc shown in two different vertical axes and the 
data are shown as a function of the silicon layer thick- 
ness. The effective refractive index will be used in the 
derivation of the impedance of the transmission lines. 

Wave-guiding impedance characterizes the wave prop- 
agation in the medium. Higher impedance corresponds 



1 



to harder movement of current carriers and lower 
impedance corresponds to easier movement of current 
carriers. The wave impedance is given by 
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where Z is the impedance, V is the voltage difference 
between metal-insulator interfaces, and / is the travel- 
ling current along the line. C is a closed contour inside 
the dielectric layer and the metal/dielectric interfaces. 
The electric and magnetic fields [Ez and -ffy), which are 
given by the solutions to the Maxwell's equations, con- 
sist of two terms due to the superposition of the surface 
waves on each metal-insulator interface. Substituting the 
resulting fields in (1) leads to 
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where ui is the operating radial frequency, W is the unit 
width of the structure, kx is the longitudinal wave vector 
of the propagating surface wave, and kzd is the trans- 
versed wave vector of the evanescent wave in the trans- 
verse plane with respect to the incident wave. For very 
thick dielectric layers (large d) or high refractive index 
material (large n), the impedance is given by 
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On the other hand, for very small dielectric layer thick- 
nesses (small d) , the impedance can be approximated by 
the first dominant term in its Taylor series, i.e.. 



^thii 



WuJtd 



(4) 



The resulting term in (4) is exactly the approximated 
formula derived previously in [8,9]. However, this ap- 
proximation by the first dominant term is only valid for 
very thin dielectric thicknesses with low refractive in- 
dex profiles. The impedance formulation as well as its 
first dominant term are shown in Fig. 2. The behavior 
of the impedance as a function of dielectric thickness is 
highly linear for very low thicknesses, where the modes 
are strongly coupled. A hand drawn sketch shows the 
behavior of the field profiles inside the dielectric layer 
in Fig. 2. Increasing the thickness of the dielectric layer 
leads to an almost constant behavior of the impedance. 
This phenomenon is due to the decoupling of the modes 
traveling on each Metal-Insulator (MI) interface. A hand 
drawn sketch shows the behavior of almost decoupled 
surface waves on MI interfaces. Fig. 2 also compares the 
proposed formulation with the other formulations pre- 
sented in the literature [8,9]. Next, we present two case 
studies for the validation of our results. 
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Fig. 2. MIM characteristic impedance as a function of the 
thickness of the silicon layer for W = Ifim. The leftmost 
shaded region is the strongly coupled region, where the 
impedance is directly proportional to the thickness of the 
dielectric layer. The middle and right shaded regions are the 
moderately coupled and the weakly coupled regions, respec- 
tively. In the weakly coupled region, the impedance is almost 
unvarying with the dielectric thickness and the modes are al- 
most decoupled. The sketch of the strongly and weakly cou- 
pled regions are also shown in the bottom. 



Tee junction design/test scenario: Tee junction 
is a power divider (shown in the inset of Fig. 3(b)), in 
which the design objective is to minimize the power that 
is refiected from the junction. The transmitted power on 
both branches of the Tee junction are equal in magni- 
tude. In order to obtain zero reflection from the Tee- 
junction, the thicknesses of the dielectric layers should 
be chosen such that the input impedance matches the 
impedance of the output branches. 

By assuming a fixed thickness of 50nm for the output 
branches, the thickness of the input layer can be calcu- 
lated using the proposed impedance formulation in (2). 
In Fig. 3(a), the reflection coefficient is shown as a func- 
tion of the ratio of the thicknesses of the input and the 
output layers. Three different methods have been used 
to derive reflection. The result from the impedance for- 
mulation in [8,9] (i.e., Zthin) is denoted by the blue solid 
line in Fig 3(a). As can be seen, our formulation, shown 
by the red dashed line, closely follows the full wave sim- 
ulation (COMSOL) results, shown by the red circles. Al- 
though the formulations in [8, 9] are relatively accurate 
for very thin dielectric layer thicknesses, the accuracy 
of the formulations deteriorates as the thickness of the 
dielectric layer increases. 

In Fig. 3(b), the optimal ratio of the input to the 
output thicknesses (i.e., where the input and output 
impedances are matched) is demonstrated as a function 
of the input branch thickness. By using the optimal ra- 
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Tabic 1. Reflection coeflicieiit of an impedance disconti- 
nuity in MIM structure 
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Fig. 3. (a) The reflection coeflicient as a function of the ratio 
of the input to the output branch thicknesses using models 
presented in the [8,9] (blue solid line), our proposed formula- 
tion (red dashed line) and the COMSOL results (red circles), 
(b) The optimized ratio as a function of the input branch 
thickness. in a Tee-junction is shown in the inset. 

tio, the Tee-junction can be designed for the minimum 
reflected power. The power flow is shown in the inset of 
Fig. 3(b) for input line thickness of 126nm and output 
line thickness of 50nm. 

Impedance discontinuity design/test scenario: 
fmpedance discontinuity is a change in the impedance 
of the line that is obtained by varying the physical size 
of the line and/or changing the dielectric material. The 
impedance discontinuity can be utilized in the design 
of Bragg reflectors. In this part, we consider the phys- 
ical thickness modulation of the line in addition to the 
dielectric material mismatch on both sides of line dis- 
continuity. Impedance discontinuity results in wave re- 
flection right at the discontinuity. The waveguide dis- 
continuity of the MIM lines is studied in [10], where the 
modal-analysis of the waveguide modes is used to analyze 
this structure using the impedance formulation in [8]. 
In the case where different dielectric layers are used at 
the sides of the discontinuity, an accurate formulation 
for the impedance is crucial for accurate calculation of 
the reflection coefficient. Based on the proposed formu- 
lation, we consider a thin MIM structure with silica as 
its insulator, which is connected to a thick MIM struc- 
ture with silicon as its insulator. In this scenario a thin 
layer of silica (e.g. 88nm thick) is terminated to a thicker 
layer of silicon (e.g. 300nm thick) both sandwiched be- 
tween silver layers. The respective impedances of the Sil- 
icon and Silica layers as well as the reflection coefficients 
are listed in Table 1. The impedance values reported by 
COMSOL are derived using the postprocessing module 
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and integrating over the desired contours. Increasing the 
number of mesh elements enhances the accuracy between 
our proposed formulation and the COMSOL results. For 
the ranges of interest, the previously reported results are 
very inaccurate and unreliable while the formulation pro- 
posed in this letter can be used to accurately obtain the 
impedance and the reflection coefficients. 

In conclusion, this letter presented a novel formula- 
tion for the characteristic impedance of the MIM trans- 
mission lines. The proposed impedance formulation is 
derived by solving Maxwell's equations. This formula- 
tion is validated using finite-element simulation software 
(COMSOL) for two test scenarios including: almost zero 
reflection from Tee-junction and the impedance discon- 
tinuity. The behavior of the impedance formulation was 
shown to be almost linear as a function of dielectric core 
thickness for very low thickness and low refractive in- 
dex materials. On the other hand, the linear behavior 
approaches a nearly constant impedance for high thick- 
ness or high refractive index materials. The almost linear 
characteristic is due to the strongly coupled modes while 
the nearly constant behavior of the impedance is due to 
the weakly coupled modes. We stress that our results 
are also applicable to moderately coupled region that 
can not be explained using prior work. 

This work was done in part when the authors were af- 
flliated with Rice University. The authors are also grate- 
ful to professor Yehia Massoud for bringing the impor- 
tance of this problem to their attention. 
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